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We start by discussing some theoretical issues of renormalization group transformations and Monte Carlo 
renormalization group technique. A method to compute the anomalous dimension is proposed and investigated. 
As an application, we find excellent values for critical exponents in \(f)z- Some technical questions regarding the 
hybrid algorithm and strong coupling expansions, used to compute the critical couplings of the canonical surface, 
are also briefiy discussed. 



1. The MCRG calculation 

1.1. Introduction 

The Monte Carlo renormalization group 
(MCRG) technique consists of the numerical 
determination of the couplings (and their deriva- 
tives) of the RG transformed action. What is 
remarkable about this method is that, if the RG 
transformation(RGT) is sufficiently short ranged, 
one may compute, in a finite lattice, critical quan- 
tities as if working on an infinite lattice. 

Rather heuristically, one may see that as fol- 
lows; suppose that a RGT possesses a Fixed 
Point(FP) consisting in a sum of local operators, 
each one just involving fields(or spins) separated 
at most by lattice sites. Assuming periodic 
boundary conditions, this FP may be accommo- 
dated in a lattice as small as (2?is -I- 1)"*. Let us 
consider it now in a (2(2ns -I- 1))'' volume, and 
apply a RGT, reducing it to a (2ns + l)*^ lat- 
tice. As the FP still fits in, one may expect the 
FP couplings (and their derivatives) not to feel 
the lattice (except possibly by tiny effects dying 
off exponentially with the linear lattice size) . Of 
course, the same line of reasoning does not fol- 
low for observables, as it is then the correlation 
length what should be fitted in, which is oo. 

It appears crucial then, both for a practical im- 
plementation and a theoretical understanding of 
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the RG, to study the properties, such as locality 
of the FP, convergence of eigenoperators, fast ap- 
proach to the FP, etc. that different RGTs have. 
This is one of the reasons for this project. 

In this work, as we are dealing with unbounded 
spins, we must determine the rescaling of the 
field, which at the FP is related to the confor- 
mal anomalous dimension rj. This is a difficult 
problem and we propose and investigate a 
variation of the Bell- Wilson criteria 

1.2. The RGT 

RGTs in real space are (exponentially)short 
ranged. The simplest of those, introduced by Bell 
and Wilson ^ is 
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There are two free parameters aw and h. 

Recall that, for unbounded spins, we may 
rescale the fields. Then, if a local FP exists, there 
is a whole line of equivalent FPs, so there is a 
marginal (eigenvalue is 1) redundant direction. 
As we stick to one of those FPs, we determine 
the parameter b by preventing moves along this 
direction. Nevertheless, when one goes to a finite 
lattice approximation, it is not clear that such a 
line of FPs still exists, the reason being that it is 
not guaranteed that all FPs are sufficiently short 



ranged. In the other case, the parameter aw just 
labels different RG transformations. It is a free 
parameter at our disposal to play with. 
Our canonical surface is 

^['^i = E^E(<^(")-'^("+A'))'+ (2) 
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and we expand the RG transformed action in 
Eq. as a sum 

i 

where the set of local operators are shown in fig. |l|. 
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Table 1 

Critical couplings at the 
bars are not written, all 
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Figure 1. Operators included in the RG trans- 
formed action. 



The first 11 coefficients of this expansion were 
tabulated in |4|. To extract critical couplings 
from those, presents some technical difficulties re- 
lated to the appearance of antiferromagnetic sin- 
gularities in loose-packed lattices, which never- 
theless, may be overcome. Critical couplings used 
are given in table 0. 

3. The algorithm 

The algorithm we used in our simulation is 
the Hybrid algorithm which has two free- 
doms, namely, the size of the leap frog step St, 
and the number of integration steps ni before the 
Metropolis test is passed. While there is not much 
room to play with St, big gains come in playing 
with 71/ . 

Indeed, in table || is reported a sample calcu- 
lation. Gradually increasing ni, autocorrelations 
reduce, with a minor penalty in CPU time. In a 
24^^ lattice, we took n; = 60 or ni = 120. 



2. The critical surface 

Following the discussion in the preceding sec- 
tion, a MCRG calculation must be done at the 
exact critical point of the infinite volume system. 
In Eq. H, for each value of A there exists a critical 
value mc(A). A strong coupling expansion, which 
for this model amounts to an expansion in the 
hopping parameter k, allows to extract all criti- 
cal couplings with negligible CPU time. 



4. Results 

We perform 3 RGTs for different aw values in a 
24^ lattice including up to 27 operators. For each 
value of aw, there exists a b such that the eigen- 
value of the T matrix in the second RGT is l(we 
discard the first RGT as we assume we are not 
yet on the linear region around the FP). If this 
value for b makes the marginal eigenvalue stable 
against the last RGT, we assume we reached a lo- 
cal FP. So, we fix b and choose the aw so that the 
canonical surface, Eq. 0, is optimally close to the 



FP. There are two cross checks. First, comput- 
ing the flow (using of Schwinger-Dyson equations 
|6|), and secondly, performing 2 RGTs in a 12^ 
lattice, and testing if expectation values of oper- 
ators in the last RGT agree with the ones coming 
from the third transformation on a 24^ lattice. 

For small values of A, though strictly we are 
in the domain of the Wilson FP, we are on the 
linear region around the Gaussian FP, and indeed, 
results are totally compatible with those. 

At the strong coupling we studied different val- 
ues of A. The most extensive one was performed 
at A = 12.672, where the optimal values where 
aw — 25 and rj = 0.03. As an example, in 
fig. H we plot the value of the second eigenvalue 
at ?7 = 0.03 for different values of aw clearly sin- 
gling out aw = 25. The couplings of the action in 
the last two RGTs agree within statistical error 
bars, though those are surprisingly large. Fur- 
thermore, the agreement between expectation of 
the third RGT and the ones that come from a 
second RGT starting in a 12^ differ, at most, by 
a 5%, which is good but not fully satisfactory yet. 

Our preliminary results for the critical expo- 
nents are 



rj = 0.030(5) 
Lo = 0.77(5) 



= 0.625(7) 
A_2 = 0.28(9) , 



(4) 



where error bars are just tentative. Those figures 
are excellent when compared with the most ac- 
cepted results rj ^ 0.035 , ~ 0.631 , oj ^ 0.8, 
(2th irrelevant eigenvalue is not known to us). 



ni 


St 


(Oi) 




CPU 


6 


0.01 


0.27449(140) 


500 


140 


6 


0.03 


0.27425( 57) 


215 


140 


6 


0.06 


0.27506( 40) 


100 


140 


11 


0.06 


0.27527( 23) 


35 


160 


21 


0.06 


0.27510( 12) 


8 


189 


31 


0.06 


0.27511( 6) 


4 


218 



Table 2 

The results correspond to a A(/)| theory in a 10"^ 
lattice at = -3.32, A = 6.0, 1.2 x 10^ configu- 
rations, throwing away 6 x 10^ for thermalization, 
Oi — 4>^ , T is the autocorrelation time, CPU is 
in units of 10"^ sec. in a Power Challenge. 
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Figure 2. Second eigenvalue dX rj ~ 0.03. dia- 
monds correspond to aw ~ oo, squares aw — 35, 
triangles aw = 25, circles aw — 20 and asterisks 
aw = 8. Dotted line is the 2nd RGT and dashed 
line the 3rd one. 
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